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Abstract. We report on a class of configurations of permanent magnets on an 
atom chip for producing ID and 2D periodic arrays of magnetic microtraps with 
non-zero potential minima and variable barrier height for trapping and manipulating 
ultracold atoms and quantum degenerate gases. We present analytical expressions 
for the relevant physical quantities and compare them with our numerical results and 
with some previous numerical calculations. In one of the configurations of permanent 
magnets, we show how it is possible by changing the angle between the crossed periodic 
arrays of magnets to go from a ID array of 2D microtraps to a 2D array of 3D 
microtraps and thus to continuously vary the barrier heights between the microtraps. 
This suggests the possibility of performing a type of 'mechanical' BEC to Mott insulator 
quantum phase transition in a magnetic lattice. We also discuss a configuration of 
magnets which could realize a two-qubit quantum gate in a magnetic lattice. 



1. Introduction 

Magnetic lattices consisting of periodic arrays of current-carrying wires [H [2] or per- 
manent magnetic films [31 H] have recently been proposed as an alternative approach 
to optical lattices O [6] for trapping and manipulating small clouds of ultracold atoms 
and quantum degenerate gases, including Bose-Einstein condensates. Magnetic lattices 
may be considered as complementary to optical lattices, in much the same way as mag- 
netic traps are complementary to optical dipole traps: they do not require (intense and 
stable) laser beams and there is no decoherence or light scattering due to spontaneous 
emission; they can produce highly stable and reproducible potential wells leading to high 
trap frequencies; and only atoms in low magnetic field seeking states are trapped, thus 
allowing the possibility of performing rf evaporative cooling in situ in the lattice and 
the study of very low temperature phenomena in a periodic lattice. Simple ID magnetic 
lattices consisting of periodic arrays of traps or waveguides have been constructed using 
both current carrying wires [3, [H] and permanent magnets P, [lOl [HI [12] on atom chips. 



In a recent paper, we proposed a 2D magnetic lattice consisting of two crossed layers of 
periodic arrays of parallel rectangular magnets In this paper, we report on a new 
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class of configurations of permanent magnets comprising four periodic arrays of square 
magnets of different thickness for producing ID and 2D arrays of magnetic microtraps 
with non-zero potential minima and variable barrier height. Analytical expressions are 
presented for various physical quantities and compared with numerical calculations. 
In one of the configurations of magnets, we show that by varying the angle between 
the crossed arrays of magnets it is possible to go from a ID array of microtraps to a 
2D array of microtraps. The 2D magnetic lattices consisting of four periodic arrays 
of square magnets of different thickness reported here may prove easier to implement 
experimentally than the crossed periodic arrays of parallel rectangular magnets proposed 
previously P]. 



2. Periodic arrays of square permanent magnets with bias magnetic fields 

Figure [Il^a)-(b) shows a configuration of four periodic arrays of square magnetic slabs 
with thicknesses ti, ^3 and ^4, respectively. In section WA\ we find the components 
of the magnetic field due to a single array of square magnetic slabs and then use our 
results to obtain the total magnetic field due to the four arrays of square magnetic slabs 
and an external bias field Bi = Bi^x + Biyij. 

2.1. Single periodic array of square magnets 

Here, we consider a single infinite periodic array of square magnetic slabs of thickness 
t, at a distance s from the plane 2 = 0, with periodicity a along the x- and y-directions 
and perpendicular magnetization [figure [T](c)]. The sum of this configuration and a 
similar one displaced by a/2 along the x-axis gives an infinite periodic array of parallel 
rectangular magnets with periodicity a along the ?/-direction and the same magnetization 
Mz and thickness t [figure[T]^d)]. For distances from the surface which are large compared 
with a/27r, the components of the magnetic field are [1] 

= (la) 

By = B^ySin{ky)e-''' {Ih) 

B, = Boycos{ky)e~''' (Ic) 

where Boy = -Bo(e^* — l)e^'', k = 27r/a and Bq = AM^ (Gaussian units). Considering the 
corresponding surface currents around the magnetic slabs in figure [I](d), we find that 
the array of blue (grey in grey style) square magnets produces the same component of 
magnetic field along the y-axis as the array of red (black in grey style) magnets. Thus, 
the y-component for the array of blue (or red) magnets can be written as half that of 
( fT6l) . Furthermore, figure [TJ^c) is symmetrical with respect to x and y, and therefore, for 
this infinite array of square magnets, the dependence of the x-component of the magnetic 
field on x is the same as the dependence of the y-component of the magnetic field on y. 
On the other hand, the 2;-component of the magnetic field should be symmetrical with 
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Figure 1. (a) and (b) Four periodic arrays of square permanent magnets with 
thicknesses ii, t2, ts and t^, respectively, (c) Array of square permanent magnets 
of thickness i, at a distance s from the plane z = 0, with periodicity a along the x- and 
y-directions and perpendicular magnetization Mz- (d) Single periodic array of parallel, 
rectangular magnets with perpendicular magnetization. 



respect to x and y and also needs to be maximum at x = y = 0. Thus the components 
of the magnetic field for the configuration of magnets in figure [U^c) may be written as 

B^ = Boism{kx)e-''^ (2 a) 

By = Bo,sm{ky)e-'' {2b) 

B^ = Boi[cos{kx) + cos(A;?/)]e"''^ (2c) 
where z :$> ^ + s + 1 and Bqi = ^{e^^ - l)e^'. 

2.2. Four periodic arrays of square magnets 

Using fl2aj) -f l2cj) . we can write the components of the magnetic field due to the four 
arrays of square permanent magnetic slabs [figure [11(a)- (b)] with external bias field 



Bi = Bi^x + Biyij as 

B^ = Box sm{kx)e-''^ + Bi^ (3a) 

By = BoySin{ky)e-''' + Biy {3b) 

Bz = [Bqx cos{kx) + Boy cos(%)]e~''^ (3c) 



Bqx — 


—Bqi + Bo2 + -Bos 


+ -Bo4 


Boy = 


Bqi — Bq2 + -Bos + 


-Bo4 


Boi = 
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where z ^ ^ + max{si + ti, S2 + ^2, S3 + ts, S4 + t^) and 

(4a) 

m 

1,2,3,4) (4c) 

The magnitude of the magnetic field above the magnetic arrays is then 
Bix,y,z) = {Bl + Bl 

+ 2[BoxBixSm{kx) + BoyBiySm{ky)]e-^' 

+ [Bl + Bl + 2Bo.Boy cos(A;x) cos(%)]e-2'=^}^ (5) 

This configuration of square magnets gives a 2D periodic lattice of magnetic traps with 
non-zero potential minima given by 

'oxBiy — -Bq, 

{Bl + Biyt^ 



p _ I Bp^Bly - BpyBlx I 



which are located at 



a^mm = (nx- + I ) a, n^: = 0, ±1, ±2, • • • (7a) 
riy + ^ ) a, = 0, ±1, ±2, • • • (76) 
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Zmin = — In [ — - — ^ ^ — ] (7c) 

27r y —Bq^Bix — BoyBiy J 

If we require a 2D lattice with non-zero magnetic field minima, ([6]) and (ITcj) impose the 
constraints 

BoxBiy 7^ BoyBix (8 a) 

Bqx + Bpy > —BqxBix — BoyBiy > (86) 

on the geometrical parameters and the components of the bias magnetic field, where 
according to fHaj) -f l4cj) . Bqx and Boy are independent of the magnetization and depend 
only on the geometrical constants ti, t2, ts, t^, si, S2, S3 and S4. 

The potential barrier heights in the three directions are given by 

4-Box,-BL, + {BoxBiy + BoyBix) 



R2 _l R2 



Bmin, j = 1,2 (9a) 



Ox "T ^Oy 

AB' = {Bl^ + Blyy^-Brmn (96) 

where Xi = x and X2 = y- Furthermore, the curvatures of the magnetic field at the 
centre of the traps can be written as 

d'^B An'^ BoxjBixj{BoxBix + BoyBiy) 



dx,^ (Bl + Biy)-2\BoxB,y - B^yB^ 



J = 1,2 (10a) 
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106) 



dz"^ dx^ dy"^ 

The trap frequencies, for an atom in hyperfine state F with magnetic quantum number 
mp, are given by 



27r7 



Bqx,Bij,^{Bq.j.Bix + BoyBiy) 



{Blx + Bly) 2 \BQxBiy - BoyB 



lx\ 



J = 1,2 



:na) 



[lib) 



where 7 = ym^gpiXg/m, is the Lande g-factor, is the Bohr magneton and m is 
the atomic mass. 



2.2.1. Symmetrical 2D magnetic lattice To create a 2D lattice that is symmetrical 
with respect to x and y we impose the constraint ^B^ = AB^ from which we obtain 
Bly = aoBix, where ao = Bo^/Boy, Bq^ 7^ and Boy 7^ 0. This is the same as the 
condition given in [5] for two crossed layers of infinite periodic arrays of magnets with 
bias fields (figure [2] (a) of [1]); so our results are the same as the equations given in [4j|. 
Here again, our analytical expressions for Bmin, Zmin, AS^, AB^, AE^ 
ujx, and ujz depend on the x-component of the bias field Bix only, rather than on 
both Blx and Biy-. 



Br, 



ai\B 



lx\ 



d^B 



a j a2Box 
— m 

271 

d^B 



\Bix\ J 
Id^B 



AB' 



UJ, 



dy"^ 

v/2 a 
ABy = aABix 



2 dz"^ 



Bu 



>lx 



AB' = a^\B 



lx\ 



where Bqx and Boy are given by (l4al) and MM and ai, 
which depend on and are given by 



a2 



\l-al\ 



2,h 



1 
2 



an 



(12a) 
(126) 

(12 c) 

(12 c/) 

(12e) 
(12/) 

as are dimensionless constants 

(13a) 
(136) 
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Figure 2. (a) Two crossed arrays of parallel rectangular magnets with perpendicular 
magnetization, (b) Three arrays of square magnets with different thickness. 
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For a symmetrical 2D magnetic lattice with non-zero potential minima above the surface 
of the top array, we have the constraints which are the same as in |1] 

a2BQx > \Bix\ > 0, aoa2Boy > \Biy\ > (14a) 
Bo. ^ Boy (146) 



2.3. Special cases of magnetic lattices 

In this section, we consider special cases of the magnetic lattices introduced in 
section I 



2.3.1. Single infinite periodic array of rectangular magnets When S3 = Si = 0, 
t^ = ti = t and ^4 = ^2 = [figure [1] (d)], we obtain an infinite periodic array of 
rectangular magnets and (l^fflj) - (l4cj) give 



5, 



Ox 



0, 



5, 



Oy 



Boie' 



kt 



(15) 



Substituting Bqx and Boy from (fTSl) into (j5l)- (111 fell , we obtain expressions for -Bmm, Zmm, 

d^B d^B d^E 



AB^, AB^, AB^, , , i^x, ^y and cuz which are the same as in P]. 



2.3.2. Two crossed arrays of parallel rectangular magnets Here, we consider a 
configuration of square magnets in which ^3 = ti, ^4 = ^2, 52 = Si = and ^4 = ^2 = s+ti 
[figure [2]^a)]. From (I4aj) -(l4cl). we have 

Box = Boie''' - l)e'^(*i+^), Boy = Bo{e^'' - 1) (16) 

Equation (fT6!) is in agreement with the expressions given in [4j| where this magnetic 
lattice has been studied in detail. 
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2. 3. 3. Three arrays of square magnets with different thickness When Si = S2 = S3 = 
^4 = [figure [2](b)], from fHaj) - fj¥cj) . we have 



We note that, according to fll4fep . for a symmetrical 2D lattice with non-zero magnetic 



potential minima, we should have t2 7^ ti- Using equations (l5|)- (]116l) and ( fT71) we 



can determine analytically the relevant quantities for this system. This magnetic 
lattice has already been studied numerically [4J . Table [1] gives the numerical ^ and 
analytical inputs and table [2] gives calculated values for the relevant quantities. There is 
good agreement between the analytical results for an infinite lattice and the numerical 
results [1] which were obtained in the central region of the lattice. 

2.3.4- Chessboard configuration of square magnets In figure [2]^b), if we take ti = t2 
and ^3 = 0, we have a chessboard configuration of square magnetic slabs. In this case, 
from ( JT71) . we obtain 

Box = Boy = (18) 



Considering fll4&p . we do not have a lattice of non-zero magnetic field minima. In fact. 



we do not have a magnetic lattice, because according to (I3aj)-(l3cj) we obtain a uniform 
magnetic field if we look at points sufficiently far from the surface of the array. 



Table 1. Numerical 4 and analytical input parameters for a configuration consisting 
of three arrays of square magnets with different thickness [figure [2] (b)] 



Parameter 


Definition 


Numerical inputs 


Analytical inputs 


Tlsq 


Number of squares 


401 


00 




in X- or y-direction 






a {^J.m) 


Period of magnetic lattice 


1.000 


1.000 


Ix = ly [tJ-m) 


Length of the chip 


200.5 


00 




along X 01 y 






ti (pm) 


Thickness of magnetic 


0.120 


0.120 




film (first) 






t2 (pm) 


Thickness of magnetic 


0.100 


0.100 




film (second) 








Thickness of magnetic 


0.220 


0.220 




film (third) 






4^M, (G) 


Magnetization along z 


3800 


3800 


Bix (G) 


Bias magnetic field along x 


-5.00 


-5.00 


Bly (G) 


Bias magnetic field along y 


-4.22 


-4.22 


Bi. (G) 


Bias magnetic field along z 


-1.87 
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(z^z'=o) 



Figure 3. Two crossed arrays of parallel rectangular magnets at an arbitrary angle 9 
with respect to each other. 



Table 2. Numerical [4] and analytical results for the input parameters from table [T] 
for a configuration consisting of three arrays of square magnets with different thickness 
[figure Kb)] 



Parameter 


Definition 


Numerical results 


Analytical results 




X co-ordinate of potential minimum 


0.250 


0.250 




Vmin iprn) 


y co-ordinate of potential minimum 


0.250 


0.250 






z co-ordinate of potential minimum 


0.952 


0.952 




d (/im) 


Distance of potential 
minimum from surface 


0.732 


0.732 






Magnetic field at potential minimum 


1.10 


1.10 




(^) 

ox ^ cm ' 


Curvature of B along x 


7.52 X 10^° 


7.52 X 


lOio 


d^B / G \ 
d'lp- ^ cm?- ' 

# (^) 


Curvature of B along y 


7.52 X 10^° 


7.52 X 


lOio 


Curvature of B along z 


1.50 X 10" 


1.50 X 


10" 


AB"" (G) 


Magnetic barrier height along x 


8.10 


8.10 




ABv (G) 


Magnetic barrier height along y 


8.09 


8.10 




AB'- (G) 


Magnetic barrier height along z 


5.45 


5.45 




Umin/kB {p.K) 


Potential energy minimum 


74 


73 




AUyks itiK) 


Potential barrier height along x 


544 


544 




AUy/ks iiiK) 


Potential barrier height along y 


544 


544 




AU'/kB ipK) 


Potential barrier height along z 


366 


366 




uJxI'i'T^ (kHz) 


Trap frequency along x 


350 


350 




LOyjlix (kHz) 


Trap frequency along y 


350 


350 




a;^/27r (kHz) 


Trap frequency along z 


494 


495 




hu>x/kB ip-K) 


Level spacing along x 


17 


17 




huy/kB i^J.K) 


Level spacing along y 


17 


17 




fiujz/kB [tJ.K) 


Level spacing along z 


24 


24 
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3. Two crossed periodic arrays of parallel rectangular magnets at an 
arbitrary angle with respect to each other 

3.1. General case 

Here, we consider a configuration of magnets which is a generalization of section I2.3.2I 
Figure [3] shows this configuration of rectangular magnetic slabs in which the upper array 
of parallel magnets is rotated at an arbitrary angle 9 about the z axis which is normal to 
the surface. The magnets in the lower array are parallel to the x-axis while the magnets 
in the upper array are parallel to the x'-axis. If 6' = 7r/2, according to figure [3l we have 
the two crossed infinite periodic arrays of magnets shown in figure [2] (a). Considering 
the symmetry, in the rotated coordinate system of reference, we have for the upper array 

B^, = (19a) 

By, = Boy, sm{ky')e-'"' (196) 

B,, = Boy' cos{ky')e-^'' (19c) 

where B^yi = i?o(e'^*^ — l)e^*^, y' = — xsin^ + ycos6 and z' = z. Now, we can write the 
components of the magnetic field in the original coordinate system 

B^ = Boy' sin 9 sm{ky')e-''' (20a) 

By = Boy, cos 9 sin(%')e"''^ (206) 

B, = Boy, cos{ky')e~^^ (20c) 

The components of the total magnetic field due to both arrays of parallel magnets plus 
a bias field Bi = Bi^x + Biyij are then 

B^ = -Boy, sin 9 sm{ky')e''''' + Bi^ (21a) 

By = [Boy/ COS 9 sm{ky') + Boy sm{ky)]e~'''' + Biy {21b) 

B, = [Boy, cos(%') + Boy cos(%)]e-'=^ (21c) 

where ^ S> ^ + ^2 + ^2 and Boy = Bo{e^^^ — 1)6^^*^ The magnitude of the magnetic field 
is then 

B{x, y, z) = {bI + Bfy + 2[Boy>{-Bi, sm9 

+ Biy COS 9) sin(%') + BoyBiy sm{ky)]e~''^ 
+ {Biy, + Biy + 2Boy'Boy[cos{ky') cos(%) 

+ siTi{ky')?AYi{ky)cos9])e-^^'Y (22) 

If > 0, these arrays of rectangular magnets can give a 2D periodic lattice of magnetic 
traps with non-zero potential minima given by 

^ _ I — BpyBix + BpyiBix COS 9 + Boy'Biy sin^ 6^1 
{Bl^, + Biy -2Boy> Boy cos 9)'^ 
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The potential minima are located at 



nx + - 



a sin 6 



A In 

27r 



Ux' + - 

1 

n,,/ 

y 4 



4/ 1-cos^' 
- I a, 



= 0,±1,±2,--- 



n„ = 0,±1,±2, ■•■ 



2BQyl B^y COS 9 



Boy'Bix sin 9 



a sin 9 



1 

a, 



cos 



Boy'Biy cos — BoyBiy ^ 
0,±1,±2,--- 



0,±1,±2, 



(24a) 
(246) 

(24c) 

(24 d) 
(24e) 



The period p of the magnetic potential minima, which are located along straight lines 
parallel to the x- and (defined in figure [3]), varies with 9 and is given by 

a sin 9 , , 

If 9 changes between and 7r/2, the period of the magnetic lattice varies between oo 
and the minimum period a, which is the period of the magnetic slabs along the y- and 
j/'-axes. Now, we can write the constraints for a lattice with non-zero magnetic field 
minima 

— BoyBix + BoyiBix COS 9 + B^yiBiy sin^ 9 Q (26 a) 



B^yi + B^y — 2BQyfBQy cos9 > 



-Boy'B 



Ix 



sin 9 



Boy'Biy cos 9 



BoyBiy > 



{26b) 



According to fl23l) - fl24cp and the contour plots of the magnetic field for different values of 
9 (figure HI), apart from ^ = for which we have a ID array of 2D magnetic microtraps, 
we have a 2D array of 3D microtraps with different periods in the x- and a;'-directions, 
as 9 varies. Looking at figure |4] more carefully, we find that as 9 changes the distance 
between the magnetic traps along the interior bisector of the angle 9 between the positive 
y- and y'-axes also changes. Thus, depending on the barrier heights, the tunnelling rate 
between the traps in this direction will be higher than the tunnelling rate in the x- and 
x'-directions. Here, we give the barrier heights in the x-, x'- and ^-directions 

-\ 2 



Bxx - (-1) 



B^yl sin 9{BQyBly — B^ylBly COS 



Boy'Bix sin( 



^Oy' 



+ 



+ ( ly ^^^y + ^^y' ^) (BpyBiy 



Bly - 2Boy>Boy COS 9 

Boy'Biy cos 9 + Boy'Bix sin 



^Oy' 



^Oy 



Bri 



2B(f yl Boy COS 9 

j = 1 or 4 



/\B' 



[Bl+BlyY^ 



Br, 



where, Xi 



and by definition we have 

^B B{x ^maxj y ymin: ^min) 



(27a) 
(276) 

(28 a) 
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-3-2-10123 -3-2-10123 -3-2-10123 



K(ym) x{|jRi) x(iJT(i) 




JC(lJM) :-Ltjjrn) x(lJin) 



Figure 4. Density plot of the magnetic field, calculated using analytical expressions, 
for two crossed arrays of parallel rectangular magnets at an arbitrary angle 9 with 
respect to each other for the input parameters from tableland (a) 6* = (b) = 7r/18 
(c) 9 = 7r/9 (d) 9 = 7r/6 (c) 9 = tt/A (f) 9 = n/'d (g) 9 = Tvr/lS (h) 9 = 47r/9 (i) 
9 = tt/2 



where 



^B""' = B{X = X'^ax^y = y'^i^, Zmin) - Bmin (286) 
AB"" = B{x^in, ymin, Z = Oo) - Bmin (28c) 



asin^^ / 3\ 

Xmax = ^ h^aj + T , n^; = 0, ±1, ±2, ■ • • (29a) 



We also have 



1-cos^ V 4 

a sin 6 
1 — cos 9 



, asin^ / 3\ , , 

^max = -7T + 7 ' = 0' ±1> ±2, ■ • • (296) 



AB"" = B{y' = ^, y = ^, Zmin) - Bmin (30a) 
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AB^' = B{y' = y = ^, Zmin) - Bmin (306) 

Furthermore, the curvatures are given by 

-BoyBl,{Bl + 7B%) cose - 2Boy>{Bl, + B',y){Bl - B%) cos2^ 
+B,yBly,{Bl - B%) cos 3^ + 2So,(3So\, + 2Bl)B,,B,ySme 
- 4Boy' {Bly, + Bly)Bi^B^y sin 29 + 2BoyBly,Bi^Biy sin 3^] (31a) 

= (j2plQ ^^y sin^ 6'[5oy (-85oj/^5iy^ + Bly,{Bl^ - Bjy) cos 61 

-5oV(SL - Si\) cos W + 45i,(So,(So\ + So,.')Si, 
+ Boy'B^^iBl - Bo/ cos 2^) sin^^)] (316) 
d^B 271^ 



where 



{BoyBly — BoylBly COS 9 + BfyylBlj; SVCi 9) [ 



{Bly{2B0y^ + Boy' cos 9{-6B^y + Boy' cos 9{-2BQy, cos 9 

+Boy{5 + COS 2^))))) + 2Boy,Bia, cos 9{BQyr — Boy COS 9) 

{-Boy + Boy' cos 9) sin 9 + 2BoyB^y,Biy sin^ 9] (31 c) 

^2 

^-72 ^ o^3^ {BoyBiy - BoyiBiy COS 9 + i?oy'-Bi3;Sin6')[ 

2B^y cos ^(-SS^^Soj// - 2So2/,3 + 2^3^ cos 9 + IBoyBoy,'^ cos ^ 

—BoyBoy/{Boy + -Boj/' cos 6') cos26') + Boyi Bi^i^^Boy^ + -Bqj/'^) 

+ Soj,Soj,'(-9 cos e + cos 3^)) sin 9] (31 d) 

52 ^ 4^2 

— = 1 {BoyBiy - Boyi Biy cos 9 + Boy>BixSm9y (31 e) 

F = \Bl, + Bly - 2BoyBoy, cos9\ (32a) 
G — \ — BoyBix + Boy'Bix COS 9 + BoyiB^y sin 9\ (326) 
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Figure 5. Magnetic field produced by two crossed arrays of parallel rectangular 
magnets at an angle = 7r/6 with respect to each other, with input parameters from 
table [3] Curves (a)-(c) show the analytical results of the magnetic field (a) along a 
line {y = ymin, z = Zmtn) parallel to the a;-axis, (b) along a line (?/' = y'^^^, z = z„„) 
parallel to the x'-axis, and (c) along a line (x — Xmin, y = ymin) parallel to the z-axis. 



3.2. Symmetrical lattice 

For 6' 7^ 0, if we have Biy = (3q{0)Bix, where I3q{6) = B^yi saiO / {B^yi cos 9 + BQy), then 
the barrier heights in the x- and x'-directions are equaL Thus, we have a symmetrical 
2D lattice of 3D microtraps in the x- and x'-directions and, as figure [5] shows, the 
microtraps in the x- and x'-directions have the same period. This figure shows our 



Table 3. Numerical and analytical input parameters for two crossed arrays of parallel 
rectangular magnets at an angle 9 = 7r/6 with respect to each other (figure [3]): (1) 
numerical inputs (2) analytical inputs 

Parameter Definition Numerical inputs Analytical inputs 



Ur Number of parallel rectangular 1001 oo 

magnets in x- or a;'-direction 

a [jim) Period of magnetic lattice 1.000 1.000 

Ix — I'x ifJ-m) Length of the magnets 1000.5 oo 

along X- or a;'-direction 

ti (fim) Thickness of magnetic 0.322 0.322 

slabs parallel to the a;-direction 

t2 (fim) Thickness of magnetic 0.083 0.083 

slabs parallel to the x'-direction 

si (fim) Distance from the plane 0.000 0.000 

z = to the lower surface of 

the magnets parallel 

to the x-direction 

S2 (/LWTi) Distance from the plane 0.422 0.422 
z = to the lower surface of 
the magnets parallel 
to the x'-direction 

AnM^ (G) Magnetization along z 3800 3800 

Bix (G) Bias magnetic field along x -4.08 -4.08 

Biy (G) Bias magnetic field along y —1.32 —1.32 

Biz (G) Bias magnetic field along z —0.69 0.00 
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analytical results for an infinite lattice, which are in excellent agreement with our 
numerical calculations in the central region of the lattice. For the numerical calculations 
we have written a Mathematica code and have used the software package Radia [T3] 
interfaced to Mathematica. Here, the quantities of interest are 



A In 

27r 



d^B d^B 



d^B 



\Bix\ 
47r^/33(g) 



\Bix\ 



Bla 



B 



lx\ 



UJr, 



2lT'~f 



(33 a) 
(335) 

(33 c) 
(33 d) 
(33 e) 



Table 4. Numerical and analytical results for the input parameters from table [3] for 
two crossed arrays of parallel rectangular magnets at an angle 9 = ti/Q with respect to 
each other (figure [21) 



Parameter 


Definition 


Numerical results 


Analytical results 




X co-ordinate of potential minimum 


0.933 


0.933 




Vmin (pm) 


y co-ordinate of potential minimum 


0.250 


0.250 






z co-ordinate of potential minimum 


1.200 


1.200 




d {fj,m) 


Distance of potential 
minimum from surface 


0.695 


0.695 




^min (^) 


Magnetic field at potential minimum 


2.68 


2.68 




f# (^) 


Curvature of B along x 


1.05 X lO^'' 


1.05 X 


IQio 


#^ (^) 


Curvature of B along y 


5.99 X 10^ 


5.98 X 


10^ 


# (^) 


Curvature of B along z 


1.65 X 10^" 


1.65 X 




AB^ (G) 


Magnetic barrier height along x 


8.32 


8.33 




AS^' (G) 


Magnetic barrier height along x' 


8.31 


8.33 




AB^ (G) 


Magnetic barrier height along z 


1.60 


1.60 






Potential energy minimum 


180 


180 




AuykB {^iK) 


Potential barrier height along x 


559 


559 




AU'^'/kB (pK) 


Potential barrier height along x' 


558 


559 




AU^ks ipK) 


Potential barrier height along z 


108 


108 




LOx/^TT (kHz) 


Trap frequency along x 


131 


131 




ujy /2Tr (kHz) 


Trap frequency along y 


99 


99 




w^/27r (kHz) 


Trap frequency along z 


164 


164 




?iujx/kB ii-i-K) 


Level spacing along x 


6 


6 




hojy/kB i^J.K) 


Level spacing along y 


5 


5 




ht^z/kB ifJ-K) 


Level spacing along z 


8 


8 
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27r7 



uo,. 



2'7r7 



Ix 



Ix 



/\B' = f3-j{e)\B,,\ 



B 



lx\ 



where 7 = ^m^gpfig/m and Pi{0), ■ 

depend on cq = B^yi / B^y and 6. Using the definitions ci = 1 — Cq, C2 
hi{6) = 1 + Co cos6' and h2{6) = 1 + Cq — 2co cos 6', we have 

Co sin 6 
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(33^) 

(33/i) 

(33«) 
(33j) 

I3^{9) are dimensionless parameters which 

' - - 1 + eg, 



Cl 



|/ll(0)|[/l2(^)]^ 

2 Co sin 6' 
2coC2 sin"^ 6 



\cih,{9)\[h2{e)]i 

cg[-8co cos6' + C2(3 + cos 26*)] sin2 6' 
|ci/ii(^)|[/i2W]i 
4co sin^ 6 



\cih,{e)\[h2m^ 

[c\ — 2cqc\ cos 9 — 4coC2 cos 29] 5 



(34a) 
(346) 
(34c) 
(34 d) 

(34e) 
(34/) 



(34^) 
(34/1) 



|/il(^)/i2(^)| 

[l + /?o(^)1^-/?i(e) 

Subject to Biy = l3o{9)Bix, the conditions for non-zero magnetic potential minima are 
Boy 7^ Boyi (35 a 



P2{9)Boy>\B,,\>0 



(356) 



According to (33/)-( l337z|) . fl34cjj) and (134 ep . for a constant Cq, it is possible to change the 
trap frequencies by varying 9 and/or the bias magnetic field Bi. Furthermore, as (133 zp . 
P3 J , p4^p and (134/tP indicate, we can adjust the magnetic field potential barriers by 
changing the angle 9, the bias magnetic field Bi, or both. 



Table [3] gives the numerical and analytical inputs and table H] gives calculated values for 
the relevant quantities. There is good agreement between the analytical results for an 
infinite lattice and the numerical results which were obtained in the central region of 
the lattice. 



A class of permanent magnetic lattices for ultracold atoms 16 



(a) (b) (c) 




-2 -1 012 -2 -1 012 01234 



x(pm) x'(ym) z(pm) 

Figure 6. Magnetic field produced by tliree arrays of parallel rectangular magnets, 
as described in section |4l The parallel magnets in the middle and the top arrays have 
an angle with respect to the parallel magnets in the lower array of = 7r/6 in (a)-(c) 
and 9 = 7r/3 in (d)-(f). The curves show the analytical magnetic field along a line 
iy = Vniin, z = Zjnm) parallel to the a:-axis, (a) and (d), along a line {y' = v'min^ 
z — Zmin) parallel to the a;'-axis, (b) and (e), and along a line {x — Xmin, y — ymin) 
parallel to the z-axis, (c) and (f). 



4. Three periodic arrays of parallel rectangular magnets 

Here, we consider a special configuration involving three periodic arrays of parallel 
rectangular magnets. The lower array (array 1) and the middle array (array 2) are at 
an angle 6 with respect to each other, similar to the crossed arrays shown in figure [3l 
while the rectangular magnets in the top array (array 3) are parallel to those of the 
middle array and remain parallel as 6 changes. The period as is twice the periods Oi 
and 02. To obtain the components of the magnetic field due to these three arrays of 
magnets plus bias field Bi = Bi^x + Biy-y, we add the magnetic field due to array 3 to 
that of arrays 1 and 2 [ fl2Tflj) -fl2Tcl)]: 

= -Boy, sin 9 sin(%')e"''^ - Boy" sin 9 sm{ky' /2)e-''^/^ + Bi^ (36a) 

By = [Boy, cos9sm{ky') + Boy sm{ky)]e~''' + [Boy„ cos9 sin{ky' /2)]e~^'l'^ + Biy (36 

B^ = [Boy' cos(%') + Boy cos(%)]e-'=^ + Boy" cos{ky' /2)e-^'''^ (36c) 

where z > f + + and Boy = 5o(e'=*i - l)e'="i, Boy' = Bq^q^^^ - and 
B^yi, = Boie^''^/'^ — This magnetic field is shown for 9 = n/6 and 9 = n/S 

in figure E] (a)-(c) and figure E] (d)-(f), respectively. The period 03 = 1 /im, which is 
twice the periods ai and 02. The thicknesses of the layers are ti = 50 nm, t2 = 40 nm 
and ts = 5 nm and the lower surfaces of the layers are at distances si = 0, S2 = 60 
nm and S3 = 105 nm from the plane z = 0. The magnetization is as before and the 
components of the bias magnetic field are Bi^ = —4.08 G, Biy = —1.39 G and i^i^ = 0. 
According to figure El we can change the separation and also the barrier height between 
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adjacent microtraps by varying 6. This might be of interest for a two-qubit quantum 
gate process [1^. 

5. Discussion and summary 

We have introduced a class of permanent magnetic lattices for ultracold atoms and 
quantum degenerate gases. The potential barriers between the adjacent microtraps can 
be altered by varying the bias magnetic field, as indicated by the relevant equations. 
Thus, these configurations of permanent magnets together with an adjustable bias mag- 
netic field could be used to study quantum tunnelling and to realize a BEC to Mott 
insulator quantum phase transition for quantum information processing in a magnetic 
lattice. 

We have shown that for a configuration of two crossed periodic arrays of parallel 
rectangular magnets with a variable angle between the two arrays, the magnetic lattice 
goes from a ID to a 2D lattice as the angle 9 between the arrays is changed from 
to 71/2. This may prove useful for loading ultracold atoms into a 2D lattice of 3D 
microtraps by initially loading the atoms into the ID lattice of 2D microtraps and then 
by changing the angle from ^ = 0to6' = 7r/2. Also, by varying the angle 9, it is 
possible to change the depth of the traps adiabatically, which suggests the possibility of 
realizing a 'mechanical' BEC to Mott insulator quantum phase transition in a magnetic 
lattice. Moreover, the period of the magnetic lattice can be adjusted as the angle 9 
is changed. This could have applications in quantum information processing and also 
in creating new artificial crystalline materials. In a more general case, it is possible to 
have two different periods and ay' = a2 in the y- and y'-directions, respectively, 

to allow for a more flexible magnetic lattice. A configuration of magnets with three 
periodic arrays of rectangular magnets having two different periods ai = 02 and a^ may 
be useful for implementing a two-qubit quantum gate [H] in a magnetic lattice. 
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